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I.  Accomplishments: 

The  intelligent  design  of  lubricants  that  perform  under  extreme  environments,  such 
as  arctic  or  space  conditions,  requires  an  understanding  of  the  how  the  tribological 
properties  are  related  to  the  molecular  structure  of  the  lubricant.  The  purpose  of  this 
grant  was  to  understand  the  relationship  between  molecular-level  perfluoropolyether 
(PFPE)  lubricant  architecture  and  its  resulting  properties  of  tribological  relevance, 
including  rheological  properties  and  chemical  stability.  We  synthesized  a  series  of 
PFPE  lubricants  varying  key  architectural  elements.  We  characterized  the  rheological 
properties  of  the  experimentally  synthesized  lubricants.  We  performed  non-equilibrium 
molecular  dynamics  (NEMD)  simulation  to  model  the  rheological  properties  and  to 
understand  the  fundamental  mechanisms  governing  the  relationship  between  molecular 
structure  and  viscosity.  We  furthermore  performed  a  multiscale  suite  of  simulations 
combining  quantum  mechanical  (QM)  calculations  of  the  reaction  rate  constants  and 
reactive  molecular  dynamics  (RMD)  to  model  the  chemical  stability  of  the  components 
as  a  function  of  architecture. 
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IV.  Summary  of  Work 

IV.A.  The  Big  Picture 


The  intelligent  design  of  lubricants  that  perform  under  extreme  environments,  such 
as  arctic  or  space  conditions,  requires  an  understanding  of  the  hovj  the  tribological 
properties  are  related  to  the  molecular  structure  of  the  lubricant.  The  goal  of  this  \work 
was  to  develop  insight  into  the  fundamental  molecular-level  mechanisms  governing  the 
relationship  between  molecular  structure  of  a  PFPE  lubricant  and  its  (i)  rheological 
properties  and  (ii)  chemical  stability.  To  accomplish  this  task,  we  performed  appropriate 
molecular  dynamics  (MD)  simulations  and  compared  them  to  an  appropriate  standard. 
For  the  rheological  properties,  we  used  non-equilibrium  molecular  dynamics  (NEMD)  to 
simulate  the  PFPE  lubricants  under  flow.  The  standard  that  we  compared  to  was  the 
experimental  characterization  of  the  lubricants  synthesized  by  Prof.  Jamie  Adcock.  For 
the  chemical  stability,  we  used  reactive  molecular  dynamics  (RMD)  to  simulate  the 
thermal  and  catalytic  decomposition  of  model  PFPE  compounds.  The  standard  that  we 
compared  to  was  quantum  mechanical  calculations.  This  process  is  schematically 
illustrated  in  Figure  1. 
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Figure  1.  A  schematic  illustrating  the  overall  methodology  used  in  this  grant. 
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IV.B.  Rheological  Properties 


We  studied  seven  PFPE  compounds  to  establish  the  relationship  between  molecular 
structure  and  viscosity.  The  structures  of  the  seven  PFPE  compounds  are  summarized 
in  Table  1.  All  PFPE  compounds,  except  compound  4,  can  be  represented  as 
variations  of  the  generic  structure 


CF3(CF2)k-[0(CF2)n]m-0(CF2)kCF3  . 

where  m  represents  the  number  of  monomeric  units  in  the  chain,  n  is  the  number  of  CF2 
units  in  each  monomer,  and  k  is  the  number  of  CF2  units  in  each  end  group.  Thus,  m 
quantifies  the  chain  length,  n  quantifies  the  monomer  size,  and  k  determines  the  size  of 
the  functional  end  group.  In  this  study,  m  assumes  values  of  2,  3,  5  and  7,  n  assumes 
values  of  1,  2  and  3,  and  k  assumes  values  of  0  and  1.  Compound  4  is  a  variation  of 
compound  2,  in  which  an  -OCF2-  monomer  has  been  inserted  into  the  middle  of  the 
chain. 

These  compounds  represent  variations  in  key  architectural  elements  including 
the  size  of  the  monomer  unit,  the  size  of  the  end  group,  and  the  chain  length. 


# 

chemical 

formula 

structure 

backbone 

atoms 

architectural 

elements 

1 

C8O4F18 

CF3-[0(CF2)2]3-0CF3 

12 

base  (m=3) 

2 

C6O3F14 

CF3-[0(CF2)2]2-0CF3 

9 

chain  length  (m=2) 

3 

C8O3F18 

CF3CF2-[0(CF2)2]2-0  CF2CF3 

11 

end  group 

4 

C7O4F16 

CF30(CF2)20CF2  0(CF2)20CF3 

11 

OCF2O 

5 

C8O3F16 

CF3-[0(CF2)3]2-0CF3 

11 

OC3F6O 

6 

C12O6F26 

CF3-[0(CF2)2]5-0CF3 

18 

chain  length 
(m=5) 

7 

C16O8F34 

CF3-[0(CF2)2]7-0CF3 

24 

chain  length 
(m=7) 

Table  1.  Summary  of  PFPE  structures  and  architectural  elements. 


IV.B.1.  The  Experimental  Standard 

In  the  Chemistry  Department  at  the  University  of  Tennessee,  Prof.  Jamie  Adcock  led 
the  team  that  experimentally  synthesized  the  PFPEs.  Viscosities  were  measured  in  the 
Chemical  Engineering  Department.  In  Table  2.,  we  report  the  experimentally  measured 
viscosities  as  a  function  of  temperature  and  architecture.  These  viscosities  formed  the 
standard  by  which  the  subsequent  NEMD  simulations  were  compared. 
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No.\T  (K) 

273 

293 

313 

333 

353 

1 

1.68 

1.16(1.18) 

0.85 

0.64  (0.66) 

0.50 

2 

0.96 

0.70  (0.67) 

0.52 

0.41  (0.43) 

- 

3 

1.56 

1.10  (1.06) 

0.80 

0.61  (0.62) 

0.47 

4 

1.34 

0.96  (0.93) 

0.71 

0.55  (0.53) 

0.43 

5 

1.69 

1.16  (1.15) 

0.84 

0.64  (0.64) 

0.49 

Table  2.  The  zero-shear-rate  viscosities  (cP)  at  various  temperatures  for  five 
compounds.  The  viscosities  from  NEMD  simulations  using  the  revised  UFF  potential  are 
in  parenthesis. 


IV.B.2.  Non-equilibrium  Molecular  Dynamics  Simulation 

In  order  to  elucidate  the  molecular-level  mechanisms  responsible  for  the  relationship 
between  PFPE  structure  and  viscosity,  we  performed  NEMD  simulations  of  the  PFPEs 
under  flow.  There  was  no  interaction  potential  in  the  literature  capable  of  correctly 
reproducing  the  experimental  viscosities.  Therefore,  we  developed  and  parameterized 
a  new  interaction  potential  which  treated  the  fluorine  atoms  differently  depending  upon 
their  local  environment.  With  this  revised  potential,  we  were  able  to  achieve  an 
excellent  agreement  between  simulation  and  experiment,  as  shown  in  Figure  2. 


Compound  No. 

Figure  2.  Comparison  of  experimental  and  simulated  zero-shear  rate  viscosities  as  a 
function  of  lubricant  structure. 
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By  measuring  the  viscosity  of  PFPE  with  different  architectural  elements  both 
experimentally  and  in  simulation,  we  can  arrive  at  several  general  statements  regarding 
the  relationship  between  viscosity  and  PFPE  architecture.  First,  longer  chain  backbone 
lengths  increase  the  viscosity.  Second,  ether  linkages  in  the  backbone  decrease  the 
viscosity.  Third,  longer  (CF2)n  units  between  ether  linkages  increase  the  viscosity. 
These  effects  can  all  be  explained  in  terms  of  chain  flexibility.  The  large  fluorine  atoms 
create  a  relatively  high-energy  barrier  for  torsional  motion  required  to  change  from  a 
gauche  to  trans  conformations.  The  presence  of  the  oxygen  significantly  reduces  this 
energy  barrier.  The  strategic  placement  of  oxygen  in  the  chain  backbone  defines  the 
size  of  the  monomeric  units  and  results  in  the  ability  to  control  of  the  viscosity  of 
compounds  that  have  the  same  chemical  formula. 

Additionally,  we  reported  the  structural  and  rheological  properties  of  four  short-chain 
PFPEs  with  common  monomeric  units  but  of  varying  chain  length  using  NEMD 
simulation  of  planar  Couette  flow.  We  found  a  power-law  index  relating  viscosity  to 
chain  length  of  2,  which  is  generally  characteristic  of  rigid  chain  molecules.  From  the 
alignment  angle,  we  observed  that  longer  chains  align  to  a  greater  degree  in  shear  flow, 
as  expected.  We  also  saw  that  longer  chains  have  a  more  pronounced  maximum  in  the 
end-to-end  distance  as  a  function  of  shear  rate,  indicating  that  longer  chains  are 
capable  of  greater  extension  (increasing  the  end-to-end  distance  at  low  shear  rates) 
and  are  more  susceptible  to  rotational  phenomena  (decreasing  the  average  end-to-end 
distance  at  high  shear  rates).  These  observations  are  commensurate  with  the 
conceptual  understanding  that  longer  chains  have  a  greater  relative  flexibility. 

Finally,  we  provided  a  quantitative  comparison  of  the  structural  and  energetic 
properties  of  relatively  rigid  PFPEs  and  relatively  flexible  alkanes  as  functions  of  chain 
length.  From  the  end-to-end  distance,  we  found  that  following  the  same  argument 
regarding  chain  flexibility,  alkanes  show  a  more  pronounced  maximum  in  the  end-to-end 
distance  as  a  function  of  shear  rate.  From  the  alignment  angle,  we  observed  that  the 
rigid  PFPEs  align  to  a  greater  extent  than  flexible  alkanes,  which  also  is  reasonable 
based  on  arguments  of  relative  chain  flexibility. 

IV.C.  Chemical  Stability 

PFPEs  degrade  according  to  two  mechanisms;  thermal  decomposition  and  catalytic 
decomposition  as  shown  in  Figure  3.  Both  reactions  involve  scission  of  the  ether 
linkage  in  the  backbone.  In  the  catalytic  decomposition,  residual  fluorine  byproduct  in 
the  lubricant  can  react  with  an  aluminum  or  iron  surface  to  form  AIF3  or  FeFa,  which 
catalytically  decomposes  the  lubricant. 

We  examined  two  model  perfluorinated  ether  compounds,  as  shown  in  Figure  4,  to 
understand  the  role  of  chemical  stability  relative  to  the  local  environment  of  the  ether 
linkage. 
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Thermal  Decomposition  CF^OCF^  -^CF^  +C/^(9 


Figure  3.  Thermal  and  catalytic  decomposition  of  perfluorinated  ethers,  color  legend; 
F-gold,  C-gray,  0-red,  Al-green. 


t  4 

CF3OCF3 

CF3CF2OCF2CF3 

Figure  4.  Model  perfluorinated  ethers,  color  legend:  F-blue,  C-gray,  0-red. 


IV.C.1.  The  Quantum  Mechanical  Standard 

In  the  absence  of  experimental  data,  we  performed  quantum  mechanical 
calculations  and  coupled  it  with  transition  state  theory  to  compute  the  reaction  rate  for 
decomposition.  We  found  a  qualitatively  different  mechanism  for  the  catalytic 
decomposition  of  the  PFPEs  than  was  given  in  the  literature.  The  old  and  new 
transition  states  are  shown  in  Figure  5. 

Using  this  new  transition  state,  we  calculate  new  reaction  rates,  as  given  in  Table  3. 
These  new  rates  show  that  at  room  temperature,  the  catalytic  decomposition  of  these 
compounds  can  be  30,000  to  100,000  times  faster  than  was  previously  thought.  Thus, 
controlling  the  chemical  resistance  to  catalytic  decomposition  is  critical. 
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Figure  5.  Left:  the  old  transition  state  in  which  the  AIF3  simply  stabilizes  the  O.  The 
new  transition  state  in  which  the  AIF3  actively  participates  in  the  reaction  as  a  fluorine 
transfer  agent. 


compound 

transition  state 
structure 

ko(s'')*10'''^ 

Ea  (kJ/mol) 

CzFeO 

old 

1.6 

228 

CzFeO 

new 

3.5 

199 

C4F10O 

old 

0.79 

197 

C4F10O 

new 

2.4 

172 

Table  3.  Reaction  rate  constants  for  the  old  and  new  transition  s 


ates. 


IV.C.2.  Reactive  Molecular  Dynamics  Simulation 

We  developed  a  coarse-grained  reactive  molecular  dynamics  (RMD)  simulation 
technique  that  allowed  us  to  match  the  quantum  mechanical  reaction  rates  for  thermal 
decomposition  in  a  large  MD  simulation.  This  new  RMD  algorithm  demonstrates  that 
we  can  correctly  reproduce  quantum  mechanical  reaction  rates  using  a  coarse-grained 
classical  RMD  algorithm.  This  will  allow  us,  in  future  work,  to  simulated  reaction  under 
flow  for  dense  systems,  which  is  not  otherwise  possible.  A  comparison  of  the 
parameterized  RMD  reaction  rate  and  the  QM  reaction  rate  is  shown  in  Figure  6. 
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Figure  6.  Comparison  of  the  reaction  rate  as  a  function  of  temperature  for  the  quantum 
mechanical  standard  and  various  parameterizations  of  the  coarse-grained  reaction 
molecular  dynamics  simulation. 

V.  Key  Findings 

The  key  findings  in  this  work  are  as  follows 

•  Five  PFPE  lubricants  were  synthesized  and  viscosities  are  reported. 

•  A  new  interaction  potential  for  PFPE  lubricants  was  developed. 

•  The  following  statements  regarding  the  relationship  between  PFPE  structure 
and  viscosity  were  validated. 

First,  longer  chain  backbone  lengths  increase  the  viscosity. 

Second,  ether  linkages  in  the  backbone  decrease  the  viscosity. 

Third,  longer  (CF2)n  units  between  ether  linkages  increase  the  viscosity. 
These  effects  can  all  be  explained  in  terms  of  chain  flexibility. 

•  A  new  transition  state  for  the  catalytic  decomposition  of  PFPEs  was 
discovered,  that  shows  that  the  decomposition  rate  can  be  100,000  times  faster 
at  room  temperature  than  what  was  previously  reported  in  the  literature. 

•  A  new  reactive  molecular  dynamics  algorithm  was  developed  that  will  allow  for 
the  reliable  and  computationally  efficient  coarse-grained  simulation  of  chemical 
reaction  (of  PFPEs  and  other  compounds)  under  equilibrium  and  nonequilibrium 
conditions. 
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